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Abstract. A boundary value problem Ve related to a third- order parabolic equa- 
tion with a small parameter e is analized. This equation models the one- 
dimensional evolution of many dissipative media as viscoelastic fluids or 
solids, viscous gases, superconducting materials, incompressible and elec- 
trically conducting fluids. Moreover, the third-order parabolic operator 
regularizes various non linear second order wave equations. In this paper, 
the hyperbolic and parabolic behaviour of the solution of Ve is estimated 
by means of slow time t = et and fast time 6 = t/ e. As consequence, a 
rigorous asymptotic approximation for the solution of Ve is established. 

partial different equations / viscoelasticity/ superconductivity/ 
boundary layer 



Diffusion et comportement ondullux dans le model linear de Voigt 

Resume. On analyse un probleme Ve des valeurs au contour relativement a une 
equation parabolique du troisieme ordre. Cette equation regie I' evolution 
unidimensionnel de beaucoup de materiels dissipatifs comme le fluides ou 
les solides visquelastiques, les gaz visqueux, les materiels superconductibles, 
les fluides incompressibles electriquement conductibles. De plus I 'operateur 
parabolique du troisieme ordre regularise divers equations non lineaires des 
ondes du deuxieme ordre. On examine dans ce travail le comportment hy- 
perbolique ou parabolique de la solution du Ve moyennant le temps lent et 
le temps rapide. En consequence, on pose une rigoureuse approximation 
asymptotique pour la solution du Ve ■ 

equations aux derivees partielles / viscollasticite / supraconduc- 
tivite 



1. Introduction 



The parabolic equation 



(1.1) 



CsU = edxxtu + (?dxxU - dttu = -f 



describes a great deal of models of applied sciences and represents a typical example of 
hyperbolic equations perturbed by viscous terms. 



According to the meaning of /, examples of dissipative phenomena related to (1.1) are: 
motions of viscoelastic fluids or solids [1, 2, 3], heat conduction at low temperature [4, 5], 
sound propagation in viscous gases [6], propagation of plane waves in perfect incompress- 
ible and electrically conducting fluids [7]. Moreover, when / = aut + sinu — j, the equation 
(1.1) is the perturbed sine- Gordon equation which models the Josephson tunnel effect in 
Superconductivity [8]. Further applications of (1.1) arise in the study of viscoelastic plates 
with memory, when the relaxation function is given by an exponential function [9]. At 
last, one remarks that even the Navier- Stokes equations for a compressible gas with small 
viscosity, in Lagrangian coordinates, can be reduced to (1.1) with / = f{ut,Ux,Uxt,Uxx) 
[10]. 

Also the meaningful analytical results concerning the qualitative analysis of (1.1) are 
very numerous and one can refer to an extensive bibliography (e.g. [10]- [18] ). In particu- 
lar, the behaviour of solutions of (1.1) when e — )■ has been analized in various applications 
of artificial viscosity method to non linear second order wave equation [19], [20]. But, from 
a physical point of view, it would be interesting to estimate the time - intervals where 
the hyperbolic or parabolic behaviour prevails, evaluating so the influence of dissipative 
causes on the wave propagation. 

These aspects are analyzed in this paper referring to the strip problem Ve for equation 

(1.1) with a linear /. The Green function G related to Vg has already been determined 
in [21] by means of a rapidly decreasing Fourier series, and its asymptotic behaviour for 
t — )■ oo has been obtained, too. 

Now, in the hypotesis of e vanishing, appropriate estimates of G by the slow time et 
and the fast time t/e will be established. As consequence, the main result is a rigorous 
approximation for the solution of the problem Vg which holds for all t < e^^' (77 > 0). 

2. Statement of the problem 

If T is a positive constant and 

D = {{x,t) :0<x <l, 0<t<T}, 
let u{x, t) the regular solution of the boundary initial value problem: 

{dxxieut + c^u) - dttu = /(x, t), (x, t) G D, 
u{x, 0) = /o(x), ut{x, 0) = /i(x), X G [0, /], 
u(0,t)=0, n(/,t) = 0, o<t<r, 

where f{x,t) is an arbitrary specified function. 

Now, denote with w{x,t) the solution of the reduced problem obtained by (2.1) with 
e = 0. To establish a rigorous asymptotic approximation for u{x, t) when e — )■ 0, we put: 

(2.2) u{x, t, e) = e"^*u;(x, t) + r(x, t, e) 
where the error r{x,t,e) must be estimated. 



By means of standard computations one verifies that r{x,t,e) is the solution of the 
problem: 

{dxxiert + cV) - dttr = F{x, t, e), (x, t) G D, 
r{x,0) = 0, rt{x,0)=0, x€ [0,1], 
r(0,t)=0, r(/,t)=0, 0<t<T, 
where the source term F{x,t,e) is: 

(2.4) F{x,t,e) = fix,t)il-e-'')+ e~''[-eXt + e\w + w^^)], 

with X = Iw + Wxx- 

The problem (2.3) has already been solved in [21] and the solution is given by: 



(2.5) r{x, t,e) = - ( di ( F{C, r, e)G(x, C, * - r) dr, 

Jo Jo 



ri ft 

di 

/o Jo 
where G{x,^, t) is: 

r, OO 

(2.6) G{x,C,t) = - Y^ Hn{t) sin7„x sin7„,f, 

n=l 

with 

—bn^t 

and 

(2.8) h = — 2 e = ge, A; = — 7n = yn. 



3. Estimates of the Green function by fast and slow times 

When e — 7- 0, two characteristic times affect the behaviour of G, i.e. t = et {slow 
time), and 6 = t/e {fast time). To point out the different contributions an appropriate 
form of G we will considered. For this, if A^ = [-^], the G- function can be given the form: 

n N OO 

(3.1) G = -{J2 + E }Hn{t) sin(7„x) sin(7„e) = Gi + G2. 

71=1 N + 1 

where, for n < N, the functions if„ are: 

-bn'^t 



(3.2) Hn{t) 



t„V(w-i -*" V(^/«)' - 1)^ 



If Q is an arbitrary constant such that: 

2d 

(3.3) 1/2 < a < 1, Na = [—], 

the term Gi can be written: 

2 A^a N 

(3.4) Giix,C,t) = - {J2 Hn{t) + J2 i^n(i)}sin(7na;) sin(7„e) • 

71=1 AfQ+1 

It is easy to prove that when 1 < n < N^ it results: 

(3.5) y^(fc/n)2 - 1 > ^ ^^^^ ; e-^"*<e-9*^ 

Otherwise, if A^^^^ + 1 < n < A^, one has N = k — (3 with < /3 < 1, and it is: 



(3.6) ^{k/nY-l> ^^f_^^~^'"" ; e-^"^* < e-2^^*/^ 



In particular, when k is an integer one has /3 = 1 and the term te~'^'^ *' "^ must be considered 
too. 

For all £ G ]0,eo] (^o < 1)) the formulae (3.5) - (3.6) allow to obtain the following 
estimate for Gy. 

(3.7) |Gi(x,e,t)| < ^0 e-"e-'?*^ + Ai e-^'^ ^-cHje^^-- ^ 

where the constants Aq A\ don't depend on e . As e < e^^^", the prevailing term in 
(3.7) is related to the slow time r = et. So, the circular component G\ is controlled by 
the slow time. On the contrary, the hyperbolic component G2 is characterized only by the 
fast time 0. In fact, let: 

^"^•^^ ^-2c/ + 7r(l-/3)' ^^-^C 

with /3 = if A; is an integer. Observing that Vn > A^ + 1, one has (see, f.i. [21]): 

(3.9) hnH{\ ± y'l - (kjnf) > cH/e, ^1 - {k/nf > e G, 
and consequently it results: 

(3.10) \G2{x,^,t)\<Gie-^ e-''"'/'. 

So, if Mo = max{Ai, Ci} and 1/2 < a < 1, the following theorem holds: 



Theorem 3.1 For all e G (0, eoK^o < 1) and (x,t) G -D, the Green function G{x,^,t) 
verifies the following estimate: 



2j- /^2a-l 



(3.11) \G{x,tt)\ < Ao e-^e""*" + Mo e-^/^e"^ */= 

where the constants Aq, Mq do not depend on e. I 

4. On the behaviour of the solution 

Now, the remainder term r{x,t,e) of (2.2) can be estimated. Referring to the function 
/ defined in (2.4), let 



(4.1) ||F|| = maxjsup |/(x, t)|,sup [|At| + e|A — ti|] }. 

D D 

Then, one has the following theorem. 

Theorem 4.1 - Let F{x,t,e) G C^{D) and let F,Fx,Ft bounded for all t. Then, the 
error term r(x, t, e) verifies the estimate: 

(4.2) \r{x,t,e)\ < k \\F\\ (e'^ tf, 
where the constants k and rj do not depend on e and r] G (0, 1/2). 

Proof- First, by (2.4)-(2.5) one deduces: 

(4.3) \r{x,t,e)\<le [ e-^"{|At(x,r)| + e|A - u|}|G(x,^,t - r)|dT 

Jo 

+1 f'\fix,T)\\l-e-'^\Gix,^,t-T)\dT. 



/o 
Further, by means of the well- known inequality [23]: 

(4.4) e"^ < [7/(ex)]^ V7>0,Vx>0 
by (3.10) and (4.3) one can deduce: 

2/Vf =-(2o— 1)7 

(4.5) M < IIFII , t^ {3/2 _4„ .•-" + (^^.),(;_^) ^ } 

So, if a and 7 are such that 3/4 < a < 1 and [2{2a — 1)]~'^ < 7 < 1, it suffices to put 

1 ^9 A/f 

(4.6) 2r] = min{{2a - 1)7 - -, 1 - a}; k = max{-Ao, / ^2wn _ J ' 



to deduce (4.2). I 

As consequence of Theorem 4.1, finally we can observe that : 

When e — 7- 0, the solution of the problem (2.1) can be approximatated by means of the 
following formula: 

(4.7) u{x, t, e) = e~^*tt)(a;, t) + r(x, t, e) 

where the error r{x,t,e) is bounded for all t < (l/e)'^. 
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